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We measured thermal conductivity k, thermoelectric power S, and dc electric conductivity a of 
Las/g.^Pr^Caa/gMnOs, showing an intricate interplay between metallic ferromagnetism (FM) and 
charge ordering (CO) instability. The change of k, S and a with temperature (T) and x agrees 
well with the effective medium theories for binary metal-insulator mixtures. This agreement clearly 
demonstrates that with the variation of T as well as x, the relative volumes of FM and CO phases 
drastically change and percolative metal-insulator transition occurs in the mixture of FM and CO 
domains. 



PACS numbers: 75.30.Vn, 72.20.Pa, 72.f5.Eb, 72.80. Tm 



Immense resurgent activities on mixed-valent man- 
ganites reveal the importance of Jahn- Teller-type elec- 
tron/lattice coupling in addition to the double exchange 
mechanism JJpQ. Another important aspect of man- 
ganites is the existence of various-scale, real-space varia- 
tion of physical properties or parameters. For example, 
even though the orthorhombicity of the average structure 
drastically decreases with the replacement of La by diva- 
lent ions in LaMnC>3, the large variation of local Mn-0 
bond lengths still remains intact Q . In high divalent-ion 
doping ranges, charge ordering results in sheet-like ar- 
rangements of Mn 3+ and Mn 4+ ions with 10-30 A length- 
scales Q. Furthermore, ferromagnetic resonance exper- 
iments have shown the presence of two types of signals 
in ferromagnetic manganites, which was interpreted as 
evidence of electronic phase separation j^]. Various ex- 
periments suggest the existence of magnetic polarons or 
mobile ferromagnetic clusters at high temperature (T), 
which could be viewed as resulting from dynamic phase 
separation [|]]|]. Recently, an electron diffraction study 
on low Curie temperature ( Tq ) manganites has revealed 
that there coexist ferromagnetic (FM)-metallic and CE- 
type charge-ordered (CO), insulating domains [|To| , |Tl[| . It 
has been emphasized that this particular type of static 
phase separation is responsible for colossal magnetoresis- 
tance in low Tq manganites. Various theoretical models 
for mixed-valent manganites also reveal the general ten- 
dency of static or dynamic electronic phase separation 
II- 

The transport properties of metal-insulator (M-I) mix- 
tures have been perennial topics for both theoretical and 
experimental condensed matter physics [jl3| . Most of the 
experiments were performed on films with deposited M-I 
mixture Jl4| or artificial bulk M-I composites prepared 
under pressure pq ]. The total electric and thermal con- 
ductivity and thermoelectric power of binary M-I mix- 
tures were successfully explained by the effective medium 
theories ||^§. 

In this letter, we report the absolute values of 



the magnetization M, thermal conductivity k, ther- 
moelectric power S, and dc electric conductivity a of 
L a 5/8-zP r :rCa3/8Mn03 with various x and T. Various 
aspects of our results are consistent with the coexistence 
of FM and CO phases, whose relative volumes change 
with both T and x, and the percolative M-I transition in 
FM-CO mixtures. Furthermore, the T and x dependence 
of a, k, and S agrees well with the (general) effective 
medium theories for M-I mixtures. 

High-quality polycrys- 
talline specimens of La5/s_ 2; Pr a ;Ca3/8Mn03 with x=0.0, 
0.1, 0.2, 0.25, 0.3, 0.35, 0.375, 0.40, 0.42, and 0.625 have 
been prepared with the standard solid state reaction. We 
fixed the Ca concentration at 3/8 because our previous 
studies showed that Tq is optimized at this particular Ca 
doping level |j) . u of all specimens with accurate geome- 
try was measured with the standard four probe method, 
and Mwas measured with a SQUID magnetometer. Both 
k and S of the representative samples (x=0.0, 0.1, 0.25, 
0.35, 0.375, 0.42, and 0.625) have been measured from 
8 to 310 K with the steady state method. A radiation 
shield was used to obtain absolute k values fl9[| . 

The systematic T-dependent and M/H curves are 
shown in Fig. 1. M/H curves were measured in H = 
2 kOe, which was carefully chosen to align FM domains 
without influencing the CO insulating phase. These re- 
sults indicate that Las/g^Pr^Cas/gMnOs (x=0.0) is, ba- 
sically, FM-metallic below 275 K, and that the ground 
state of Pr 5 / 8 Ca3/ 8 Mn03 (:r=0.625) is CO-insulating be- 
low 225 K. The behaviors of a and M/H for other x com- 
positions are systematically in-between those for ;c=0.0 
and 0.625. Open circles (Fig. 1 (a)) represent the M- 
I transition points where a becomes (local) minimum 
p0| . Open circles in Fig. 1 (b) show the M/H values 
at the M-I transition points determined from the curves. 
The average of those M/H values (dotted line) is 17±2 
% of 8.1 emu/mol, the saturated M/H value of :e=0.0. 
Thus, with changing T, the M-I transition occurs when 
M of each sample becomes about 17 % of that of x=0.0, 
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independent from x value. If we assume that the T- 
dependent volume fraction f(T) of the FM domain is 
proportional to M(T) in H=2 kOe, the M-I transition 
with changing T occurs when / reaches ^0.17, close to 
the three-dimensional percolation threshold (/ c ). This 
behavior can be also seen in the variation of a with 
x for fixed T. In Fig. 2 (a), a vs. M X (T)/ M .o(T) 
plot is shown at 10 and 100 K for various x. Inter- 
estingly, the samples with M X (T)/M , (T) < 0.15 are 
in the insulating regions of the curves in Fig. 1, while 
those with M X (T)/ M .o(T) > 0.17 are in the metallic 
regimes. Therefore, with variation of x, the M-I transi- 
tion occurs when M x /Mq,q ssO.15-0.17. These observa- 
tion clearly suggest that the M-I transition with both T 
and x takes place when f(T,x) (=M X (T)/M . (T)) be- 
comes close to f c . It is important to note that Mqq(T) 
contains the natural T-dependence of ferromagnetic mo- 
ment so that f(T,x) equals to M X {T)/M (T) (not 
M x (T)/M o , o (T=0)). 

To gain further insights into the nature of the M-I tran- 
sition, we measured T-dcpcndent k and 5 as shown in 
Fig. 3. First, we note that the estimated K e from ct, using 
the Wiedemann- Franz law, is by two orders of magnitude 
smaller than the measured k( T) for all x, indicating the 
dominant phonon contribution. Furthermore, at high T 
above T c or T co , k always increases when T is raised, 
and the magnitude of k is in the range of 0.5-2 W/mK, 
comparable to that of amorphous solids plj . This be- 
havior has been attributed to local anharmonic lattice 
distortions associated with small polarons [^2|. Related 
to this, S(T) of all the samples at high T follows the 
form So+Eg/ksT with the gap energy E g systematically 
increasing from 4 meV (i£=0.0) to 12 meV (x=0.625). 
These E g values are significantly smaller than the acti- 
vation energies (125 meV: x=Q.Q to 175 meV: i=0.625) 
associated with ct( T). This difference can result from the 
small polaronic transport |23f| . 

As shown in Fig. 3 (a), with increasing x, k(T) 
smoothly evolves from that of x=0.0 to that of £=0.625, 
and the k increase at Tq becomes smaller. Consistently, 
k vs. M x (T)/M ,o(T) at 10 and 100 K (Fig. 2 (b)) 
shows that k varies monotonically from the maximum (k 
of £=0.0) to the minimum (k of :r=0.625). For x=0.0, k at 
Tc increases suddenly, probably due to the suppression 
of local lattice distortions associated with small polarons. 
For x=0.625, k(T) behaves as in amorphous solids in the 
entire T range, except for a slight increase at T co . 

In comparison with k and a, S exhibits seemingly dif- 
ferent behaviors with x. S is very close to the metallic 
value, S of x=0.0, even near f c where a (or n) is still sig- 
nificantly smaller than a (or k) of x=0.0. S vs. M x / Mq.o 
at 100 K (Fig. 2 (c)) clearly demonstrates this tendency; 
S is close to zero (slightly negative), and is insensitive 
to M x /Mq,q as long as M x /Mq,q ~10 %. In contrast, a 
for M x /M .o = 10 % at 100 K is more than three or- 
ders of magnitude smaller than that of :r=0.0. In fact, 



T-dependence of S near Tq is also consistent with this 
metallic S behavior near f c . With decreasing T near Tc, 
S starts to decrease, i. e., becomes metallic at T higher 
than those for k and ct changes. For example, in the heat- 
ing curves of x=0.35, S starts to decrease around 130 K, 
significantly higher than T (100 K) for abrupt k increase 
or T (~110 K) for a minimum. 

When thermal/electronic transport properties of our 
systems are viewed as those of M-I mixtures, the above 
peculiar S behavior is, in fact, consistent with the theo- 
retical prediction of effective thermoelectric power Se by 
Bergmann and Levy |l8|. For an isotropic binary mix- 
ture, they showed that in terms of a, k, and S of each 
component, Se is given by 

5e = <Sm + (Si — Sm) ( E { M — 1 J / ( — 7 M — 1 ) , 

(1) 

where the subscripts M and I refer to metallic and in- 
sulating components, respectively, ke and cte refer to ef- 
fective thermal and electric conductivity, respectively, of 
the binary mixture. This equation has been successfully 
applied to explain S behaviors of binary Al-Ge films [ p"4[ . 
When (Tj/ctm << k\ /% <1 (which applies to our sys- 
tem) and for /= f c , the above equation leads to Se ~ Sm, 
which explains our experimental results noted above. 

To quantitatively compare our results with Eq. (1), 
we calculated Se at every T In this comparison, 

experimental a and k, shown in Figs. 1 and 3, were 
used for oe and ke , and (at, «i, and Sj) and (ctm, km, 
and 5m) arc assumed to be identical with those of ^=0.0 
and 0.625, respectively. (For large x, it was difficult to 
measure 5 at low T due to high resistivity, so we assumed 
that Si changes as 1/ T, and that o\ does exponentially.) 
The solid lines in Fig. 3 (b) depict the calculated Se, 
using Eq. (1), for heating curves of :r=0.25, 0.35, and 
0.42. The calculated curves match with our experimental 
5 surprisingly well at all T below Tc or T co . 

For cte (or ke) of a binary M-I mixture, Mclachlan |ll| 
proposed the general effective medium (GEM) equation, 

\ct/ +Aa^ ) \°m + A(T e ) 

where A =(l-/ c ) // c . The same equation also works for k. 
The critical exponent t is close to 2 in three dimension. 
This equation has been successfully applied to isotropic 
inhomogeneous media in wide / regions including perco- 
lation regime 0^,0. 

To apply the GEM equation to ct(T), we assumed that 
f(T,x) = M X (T)/M . (T), ct m (T) = a(T) of x=0.0, and 
ai(T) = a(T) of ^=0.625. With the parameters t=2 & 
/ c =0.17, the calculated cte for various x are shown as solid 
lines in Fig. 4. At T > ^80 K, ct e (T) nicely matches 
the experimental cr(T) even if a changes by 6 orders of 
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magnitude with T and x. However, this agreement does 
not hold at very low T. The calculated 0e(T) at T < 
80 K with the same parameters t—2 & / c =0.17 signifi- 
cantly deviated from the experimental cr(T). We found 
that at T < 80 K, the calculated cte(T) matches the ex- 
perimental <j( T) better when t is increased to ~ 4. This 
is more evident in the x dependence of a at 10 and 100 K, 
as shown in Fig. 2 (a). cte(100 K), calculated with t=2 
& / c =0.17 (solid line), matches the experimental values 
(open circles) better than that with t=4 & / c =0.17 (dot- 
ted line). However, cte(10 K), calculated with i=4 & 
/ c =0.15 (solid line), is closer to the experimental a (solid 
circles) than that with i=2 & / c =0.15 (dashed line) [fj5| . 
(The change of / c in the range of 0.15-0.17 makes lit- 
tle difference.) These observations demonstrate that t, 
normally close to the three dimensional exponent of 2, 
becomes ~4 at very low T. A similar, drastic increase 
of t has been noted in the case of tunneling transport 
for M-I mixtures |l4"|]2lj| ], suggesting that the tunneling 
process between FM domains is important for a of our 
system at very low T [p7[ . 

By using the GEM equation for k~e, the k(T) for 
various x can be calculated with the assumption that 
ki(T)=k(T) of :r=0.625, k m (T)=k(T) of x=0.0, t=2, 
and / c =0.17. The solid lines in Fig. 3 (a) represent the 
estimated ke for x=0.1, 0.25, 0.35, and 0.42. In addition, 
the calculated k~e as a function of M x /Mq.q at 10 and 100 
K is depicted as solid lines in Fig. 2 (b). Estimated ke 
lines in Figs. 2 and 3 coincide with the experimental data 
well [ p8| . In order to confirm self-consistency, Se at 10 
and 100 K is evaluated by using the calculated cte and ke 
(solid lines of Fig. 2 (a) and (b)), and the Eq. (1). The 
calculated Se (solid lines of Fig. 2 (c)) with the variation 
of x is in good agreement with the experimental values. 

The unambiguous agreement between the measured 
thermal/electronic transport properties and the calcu- 
lated values based on Eqs. (1) and (2) strongly indicates 
that: (1) transport properties are dominated by ther- 
mal/electrical conduction in M-I mixtures, (2) the rel- 
ative volume of the (FM) metallic phase is proportional 
to the measured M(T,x), and (3) the T-dependent trans- 
port and magnetic properties of metallic and insulating 
phases are always that of x=0 and 5/8, respectively. 
Combined with the earlier electron diffraction results 
, this successful agreement demonstrates that all of 
the thermal, electronic, and magneto-transport proper- 
ties of Las/g-zPr^Caa/gMnOa are dominated by the per- 
colative conduction through FM-metallic domains which 
is statically mixed with CO insulating domains. One sur- 
prising indication from our results is that at least in low 
Tc materials (x > 0.25), the so-called Curie transition is, 
in fact, the M-I transition across percolation threshold, 
and the ordered FM moment changes smoothly near the 
percolative phase transition T. 
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FIG. 1. (a) a(T) of Las/g^Pr^Cag/gMnOs for cooling 
(dotted lines) and heating (solid lines), (b) M/H curves 
with zero field cooling. Open circles depict the M-I transi- 
tion points determined from the data. T-dependent volume 
fraction of the FM domains, f(T), for each x, was determined 
as/(7» = M X (T)/M . (T). 



FIG. 2. S, k, and a values vs. M x /M .o at 10 and 100 
K. The lines depict the theoretical predictions by Eqs. (1) 
and (2). The solid lines in s show the theoretical results with 
t=4 & / c =0.15 at 10 K and t=2 & / c =0.17 at 100 K. The 
dotted and dashed lines represent the theoretical results with 
t=2 & / c =0.15 at 10 K and with t=4 & / c =0.17 at 100 K, 
respectively. For the theoretical predictions (solid lines) of k 
at 10 and 100 K, t=2 & / c =0.17 were used. 



FIG. 3. (a) k(T) for cooling (crosses) and heating (solid 
circles). The solid lines show the predictions of n(T) 
(heating) by Eq. (2) with t=2 & / c =0.17. (b) S(T) of 
Las/g.j.Pr^Cas/gMnOs. Crosses represent cooling curves for 
2=0.42 and 0.35, and the other symbols represent heating 
data. (S of x=0. 1 is very close to that of x=0. and 0.25 
below Tc and omitted for clarity.) The solid lines show the 
predictions of Eq. (1) for £=0.25, 0.35, and 0.42 with heating. 



FIG. 4. Theoretical predictions (solid lines) of u(T) by the 
GEM equation (Eq. (2)) with t=2 & / c =0.17 above -80 K. 
Solid squares are the experimental data redrawn from Fig. 1. 
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